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Gravity Probe B, launched 20 April 2004, is a space experiment testing two fundamental predic- 
tions of Einstein's theory of General Relativity (GR), the geodetic and frame-dragging effects, by 
means of cryogenic gyroscopes in Earth orbit. Data collection started 28 August 2004 and ended 
14 August 2005. Analysis of the data from all four gyroscopes results in a geodetic drift rate of 
—6, 601.8 ± 18.3 mas/yr and a frame-dragging drift rate of —37.2 ±7.2 mas/yr, to be compared with 
the GR predictions of —6,606.1 mas/yr and —39.2 mas/yr, respectively ('mas' is milliarc-second; 
1 mas = 4.848 x 10~ 9 rad). 



INTRODUCTION 

In 1960, L. I. Schiff 1 showed that an ideal gyroscope 
in orbit around the Earth would undergo two relativistic 
precessions with respect to a distant inertial frame: 1) a 
geodetic drift in the orbit plane due to motion through 
the space-time curved by the Earth's mass; 2) a frame- 
dragging due to the Earth's rotation. The geodetic term 
matches the curvature precession of the Earth-Moon sys- 
tem around the Sun given by W. de Sitter in 1916 [2]. 
The Schiff frame-dragging is related to the dragging of 
the orbit plane of a satellite around a rotating planet 
computed by J. Lense and H. Thirring in 1918 [3]. Frame 
dragging has important implications for astrophysics; it 
has been invoked as a mechanism to drive relativistic jets 
emanating from galactic nuclei |4j. 

The measurement requires one or more gyroscopes ref- 
erenced to a remote star by an onboard telescope. In the 
642 km polar orbit of Gravity Probe B, the two effects 




-6,606 mas/year NS 



are at right angles, as in Fig. [T] The predicted geodetic 
drift rate is —6,606.1 mas/yr; the frame-dragging drift 
rate with the chosen star IM Pegasi is —39.2 mas/yr. 

GP-B was conceived as a controlled physics experiment 
having mas/yr stability (10 6 times better than the best 
modeled navigation gyroscopes) with numerous built-in 
checks and methods of treating systematics. Three prin- 
ciples guided the design: 1) make Newtonian gyro drifts 
•C the predicted GR effects; 2) add sensors so that mod- 
eling hinges on physical understanding as against the 
ad hoc observational modeling used in navigation gyro- 
scopes; 3) exploit natural effects such as stellar aberra- 
tion in calibrating the instrument. Meeting the many 
mechanical, optical, and electrical requirements rested on 
a conjunction of two technologies, space and cryogenics. 

Operation in space separates the two effects, increases 
the geodetic effect 12.4 times as compared to a gyroscope 
at the equator, eliminates 'seeing' in the measurement 
to the guide star, and vastly reduces torques from sus- 
pending the gyroscope against 1 g gravity. Cryogenics 
brings new levels of magnetic shielding, thermal isola- 
tion, ultra-high vacuum operation, and a uniquely ef- 
fective gyro readout based on the London moment in a 
spinning superconductor. The two together give the in- 
strument ultimate mechanical stability: in zero g, there 
is no sag; at zero K, there is no thermal distortion. 

Essential to GP-B as a controlled physics experiment 
was the calibration phase, a 46-day period following the 
main science phase designed to set limits on a range of 
potential disturbing effects and quantify any that might 
prove larger than expected. 

EXPERIMENT DESCRIPTION 



FIG. 1: Predicted drift rates of GP-B gyroscopes. 

See [5] for definitions of WE and NS inertial directions. 



The heart of the instrument was a 0.92 m long fused 
silica structure containing four gyroscopes and a star- 
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TABLE I: The seven near zeroes 
Property Requirement Achieved 

Rotor Properties 



Mass unbalance, Sr/r 


6 x 10" 7 


2 - 5 x 10 


Asphericity (nm) 


55 


< 33 


Patch dipole (C-m) 


< 10~ 15 


see text 


Environment 






Cross-track acceleration (g) 


lO" 11 


lO" 11 


Gas pressure (torr) 


HT 11 


< lO" 14 


Rotor trapped field (fj,G) 


9 


0.2-3 


Mixed 






Rotor electric charge (electrons) 


10 s 


< 10 8 



tracking telescope mounted in a 2440 £ superfluid helium 
dewar operating at 1.8 K. Each gyroscope comprised a 
38 mm diameter niobium-coated fused quartz sphere sus- 
pended electrically, spun up by helium gas, and read out 
magnetically. The four were set in line, two spinning 
clockwise and two counterclockwise, with axes initially 
aligned to the boresight of the telescope. The space ve- 
hicle rolled with 77.5 s period about the line of sight to 
IM Pegasi, which was, however, occulted by the Earth for 
almost half each orbit (Fig. [IJ. During occulted periods, 
pointing was referenced to star trackers and rate gyros on 
the outside of the spacecraft. Drag compensation, origi- 
nated by G. E. Pugh in an independent proposal for an 
orbiting gyroscope experiment [BJ, two months prior to 
Schiff 's paper, was by a control system referred to one of 
the gyroscopes as a proof mass. Attitude, translational, 
and roll control authority was provided by the helium 
boil-off gas from the dewar vented through proportional 
thrusters [7]. 

Vital was a gyro readout that did not disturb the spin 
orientation. Superconductivity supplied three essentials. 
The spinning rotor generated a London moment equiva- 
lent at 80 Hz to a uniform 5 x 10 _5 G field aligned with 
the spin axis. A SQUID magnetometer coupled to a su- 
perconducting loop on the gyro housing provided a read- 
out capable of detecting a 1 mas change of spin direction 
in 10 hr. Finally, a combination of high permeability and 
ultra-low field superconducting shields around the instru- 
ment with local shields for each gyro achieved: 1) 240 dB 
isolation from external magnetic disturbances; 2) a limit 
on trapped fields in the rotors ~ 1% of the London mo- 
ment; and 3) virtual elimination of magnetic torques. 

The gyro readout scale factor C g (t) was calibrated on 
orbit against the aberration of starlight. During each 
half-orbit when IM Pegasi was visible, the telescope 
pointed at its apparent position, and the 20.49586 arc-s 
annual and 5.18560 arc-s orbital aberrations, derived re- 
spectively from the JPL Earth ephemeris and GPS orbit 
data, appeared in the gyro readout. 

An elementary calculation captures the stringencies of 
the experiment. Consider a spherical, not quite homo- 



geneous rotor under transverse acceleration /. Let Sr 
separate the centers of mass and support, v s ~ 9 m/s 
be the rotor's peripheral velocity, f2o the maximum al- 
lowed drift rate 0.1 mas/yr (1.5 x 10 _17 rad/s). Then, we 
have Sr/r < 2v s flo/5f, and mass unbalance requirements 
Sr/r: 6 x 10~ 18 on Earth; 6 x 10~ 10 in a typical 642 km 
altitude satellite; and 6 x 10~ 7 for GP-B in 10 -11 g drag- 
free mode. Without drag compensation, the experiment 
would have been impossible. Likewise for an aspherical 
homogeneous rotor, the electrical support torques, while 
greatly reduced on orbit, only reached the desired level 
through a further symmetrizing factor, spacecraft roll. 
The net allowed asphericity was 55 nm. 

Aiming for a ~ 0.5 mas/yr mission, we created an error 
tree setting limits on 133 disturbing terms, and verified in 
advance that all were negligible. Central were seven 'near 
zeroes' (Table [j}. Four noted already are inhomogeneity 
and asphericity for the rotor, residual acceleration and 
magnetic field for the environment. The others are gyro 
electric charge, gas pressure, and patch effect. Six met 
requirements; the complication discovered during the cal- 
ibration phase was the patch effect. Even so, the idea of 
seeing relativity in the 'raw' data was preserved. Fig. [2] 
shows the NS drifts of the four gyroscopes with no torque 
modeling. The geodetic drift is visible in all four. 

The term 'patch effect' [8j [9] refers to contact poten- 
tial differences between crystalline or contamination re- 
gions on a metal surface. Prelaunch studies focused on 
eddy current losses and torques from interaction of the ro- 
tor's patch-induced electric dipole moment with support 
voltages and the housing dipole moment. Far more im- 
portant was the on-orbit discovery of forces and torques 
caused by interacting patches on the rotor and housing. 
Put simply, rotor and housing were spherical mechani- 
cally; electrically, they were not. 

Three unforeseen effects emerged, differing in detail for 
each gyroscope: a changing polhode path and two patch 
effect torques. The changing polhode, originating in a 
< 1 pW dissipation of the rotor's kinetic energy, compli- 
cated the C g determination. Beginning with their body 
axes arbitrarily oriented and a 2 - 6 hr period, T p (t), the 
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FIG. 2: North-South gyro orientation histories with no mod- 
eling of torque or scale factor 
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rotors transitioned to a final 1 - 4 hr period, each spinning 
nearly about its maximum inertia axis. The two torques 
were 1) a spin-to-roll misalignment torque 200 — 500 x 
larger than predicted from mechanical asphericity, and 
2) a 'roll-polhode resonance' torque, where the roll av- 
eraging mentioned above would temporarily fail, making 
a particular gyroscope axis realign, or step over, by as 
much as 100 mas in 1-2 days when a high harmonic of 
its polhode rate came into resonance with spacecraft roll. 
The patch effect misalignment torque was discovered and 
quantified during the calibration phase by commanding 
the spacecraft to point to a series of positions at known 
angles to the guide star. 

Remarkably, a key to modeling all three effects was 
magnetic asphericity: two patterns, magnetic fluxons and 
voltage patches, remained locked together in the rotor. A 
process of Trapped Flux Mapping (TFM) allowed exact 
tracking of the evolving polhode phase and angle needed 
for computing both C g (t) and torques. 

Further extensive and diversified post-science calibra- 
tions showed no significant disturbing effects other than 
the three just discussed. 



DATA ANALYSIS 

Two data analysis methods were used to determine 
and cross-check the relativity results. One, called 'al- 
gebraic', was based on a parameter estimator utilizing 
the gyro dynamics and measurement models detailed be- 
low. It produced the primary science results (see Ta- 
bles [nj[inj and Fig. [3b. The other method, called 'geo- 
metric', though currently less precise, neatly eliminated 
the need to model the misalignment torque. This torque, 
being Newtonian, causes drift rates perpendicular to the 
misalignment vector fl. The drift in that direction is 
thus a combination of relativistic and Newtonian contri- 
butions, while the drift rate parallel to fl - when there is 
no roll-polhode resonance - is pure relativity. The annual 
variation in misalignment direction allows determination 
of both relativity effects. 



third term of Eq. (JlJ where A(j> m (t) = cf> r (t) —m<f> p (t), cf> r 
and <j) p being the known roll and polhode phases, with 
exact resonance occurring when the corresponding fre- 
quencies coincide, uj t = muj p . To clarify the picture of a 
resonance, we keep only the resonance term on the right 
of Eq. 0. Then the 'step-over' in the NS— WE plane 
follows, to lowest order, a Cornu spiral winding out from 
its initial direction a day or so before the resonance, mov- 
ing across, then winding back in to the new direction, its 
angular magnitude proportional to yj a^ n + [10] . 

Science Data. Determining r-^s arL d r^E requires: 
1) data from gyroscope, telescope, and roll reference 
signals (preprocessed and synchronized at 2 s inter- 
vals), with supporting TFM analyses based on separate 
2200 sample/s gyro readout signals; 2) the measurement 
model of Eq. ^ below. The data came from 11 months 
of science operations, where spacecraft anomalies from 
events such as solar-flare-induced CPU reboots resulted 
in ten distinct but connectable data segments. Minor in- 
terruptions came from calibrations, transient electronics 
noise, and passages through the South Atlantic Anomaly. 

Measurement Model. With s, f, <p r and C g (t) mean- 
ings already defined, we write SQUID signal z(t) as: 



Gyro Dynamics. The r^g (geodetic) and r\vE (frame- 
dragging) drift rates and two patch effect torques [TU] 
connect to the unit vector s(t) along the gyro spin axis 
by the equations of motion: 

SNS = HsrS + k(t)fl WE + ^ (° m C0S ^ m _ bm sm ^4>m) 

m 

%E = f W E - fc(*Ws + X! (° m Sin Z ^ m + bm C0S A0m)(l) 

m 

The fc(i)/i\vE, k(t)fj,yis terms give the misalignment drift 
rates, /xwe, Mns being the components of the misalign- 
ment vector fx = t — s, and k(t) the polhode-dependent 
torque coefficient; f is the unit vector along the space- 
craft roll axis. The resonance torques are governed by the 



z(t) = C g (t) [(r NS - s NS ) cos(0 r - 
(twe — %e) sin((j> r + 5(f>)} + bias - 



noise; 



(2) 



(tns — s ns) and (tew — sew) are the NS and EW projec- 
tions of the misalignment ft = t — s, and 6(j> is a constant 
roll phase offset. 

Trapped Flux Mapping. Central to the three main 
modeling challenges (determining C g (t) and handling the 
misalignment and resonance torques) was the informa- 
tion gained from Trapped Flux Mapping. The magnetic 
flux trapped in each rotor on cooling through its transi- 
tion temperature formed a unique, stable fluxon pattern. 
By fitting to the spin harmonics of the 2200samplc/s 
data stream, we constructed detailed maps from which 
the three-dimensional orientation of each gyroscope could 
be tracked continuously over billions of turns. Rotor spin 
speed and spin-down rate were determined to 10 nHz and 
lpHz/sec, respectively QTj. Crucial for data analysis 
were: 1) polhode phase 4> p (t) and angle 7 P (i) both good 
to 1° over the 353 days of science; and 2) detailed knowl- 
edge of the trapped flux contribution to C g (t) . 

TFM and C g (t). The scale factor C g (t) to be calibrated 
against aberration comes from the London moment Ml 
aligned with the rotor's spin axis, plus a complex pat- 
tern of trapped flux ~ 10~ 2 Ml fixed in its body, i.e. 
C g M + C g rF (t). To meet the necessary < 10~ 4 accuracy 
for C g (t) required connecting data from ~ 5200 succes- 
sive guide star valid half-orbits enabled by TFM. 

TFM and Resonance Torques. The role of TFM in 
computing the resonance torque was even more remark- 
able. Solving Eq. fl]) required precise knowledge of 4> p (t) 
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and 7 P (i). It is curious to reflect that if GP-B had had 
an ideal London moment with no trapped flux, the com- 
putation would have been impossible. 

Spacecraft Pointing. To model the misalignment 
torque, the accurate value of f (t) is needed continuosly. 
During the half-orbit periods when the guide star is vis- 
ible it is the sum of annual and orbital aberrations and 
minor terms including parallex, bending of stra light by 
the Sun, and pointing error from the telescope signal. 
For the occulted periods, r(t) was determined from the 
four science gyros' SQUID signals. The connection pro- 
cess required fast computation of the gyro spin vector 
s(t) through the entire 11 months of science data. 

Parameter Estimator. The science data set was 
processed by a nonlinear batch weighted least-squares 
(WLS) estimator implemented in iterative sequential in- 
formation filter form [12] . The cumulative information 
matrix for the full mission (~ 1.4 x 10 7 data points per 
gyroscope) was calculated sequentially based on one-orbit 
data batches. To ensure robust convergence of numerous 
runs with different numbers of estimated parameters, wc 
used the sigma-point filter technique [13]. The lack of 
convergence observed in some runs of a standard WLS 
estimator was thus eliminated. Computations were crit- 
ically facilitated by the replacement of differential equa- 
tions with their analytical solution sns(0j s we(£) for 
an arbitrary k(t). 



RESULTS AND CONCLUSIONS 

The four gyroscope signals were analyzed indepen- 
dently and the results combined and cross-checked in var- 
ious ways. Table [n] lists the four relativity estimates and 
final joint result with 1 a errors, also plotted as 95% con- 
fidence ellipses in Fig. [3j The values are referenced to 
inertial space by correcting for guide star proper motion, 
27.3 ± 0.1 mas/yr NS and -20.0 ± 0.1 mas/yr WE [H], 
and the solar geodetic contribution of 7.3 ± 0.3 mas/yr 
NS and -16.2 ± 0.6 mas/yr WE. 

The result for each gyroscope is the average of 10 drift 
rate estimates based on 10 distinct parameter sets deter- 
mined by the following uniform procedure. The submod- 
els for scale factor C g (t), misalignment torque coefficient 



Source 

Gyroscope 1 
Gyroscope 2 
Gyroscope 3 
Gyroscope 4 
Joint (see text) 



TABLE II: Results 
r NS (mas/yr) 

-6, 588.6 ±31.7 
-6, 707.0 ±64.1 
-6, 610.5 ±43.2 
-6, 588.7 ±33.2 
-6, 601.8 ± 18.3 



rwE (mas/yr) 

-41.3 ±24.6 
-16.1 ±29.7 
-25.0 ± 12.1 
-49.3 ± 11.4 
37.2 ± 7.2 
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FIG. 3: North-South and West-East relativistic drift rate 
estimates (95% confidence) for the four individual gyroscopes 
(colored ellipses) and the joint result (black ellipse). 



k(t), polhode phase <p p (t), etc., are linear combinations of 
certain basis functions with coefficients to be estimated. 
The number of terms in each submodel is increased from 
zero until the change in the relativistic drift rate esti- 
mates becomes less than 0.5 a. This defines the baseline 
set for a given gyroscope. Table [TT] gives the weighted 
mean of the baseline and nine 'sensitivity run' estimates 
for each gyroscope, whereby in each run the number of 
terms in a single submodel is increased by 1 above the 
baseline value. The post-fit residuals were white with 
X 2 ~ 1 for all gyroscopes. 

The joint 4-gyro result is a weighted average of the 
four individual drift rates using their 2x2 tns > r WE sub- 
blocks of the full covariance matrix for each gyroscope. 
In the NS direction the scatter of the individual estimates 
is 28% larger than the individual uncertainties indicate, 
while in the WE direction it is 37% smaller. The individ- 
ual and combined statistical uncertainties are corrected 
for their 'over' and 'under' dispersion using the \ 2 °f the 
individual estimates in the NS and WE directions. The 
'parameter sensitivity' covariance matrix is the scatter 
of all 10 4 combinations of drift rate estimates obtained 
in 10 sensitivity runs for each of the 4 gyroscopes. The 
total covariance matrix is a sum of three matrices: the 
corrected statistical one, the sensitivity one, and the co- 
variance matrix of unmodeled systematic effects. The 1 a 
uncertainties of the joint result shown in Table II are the 
square roots of the diagonal elements of the total covari- 
ance matrix. 



GR prediction 



-6,606.1 



-39.2 



Table III summarizes uncertainties of the joint 4-gyro 
result. The three main disturbances (evolving scale fac- 
tor and the two Newtonian torques) are modeled and thus 
contribute to the statistical uncertainty. The systematic 



uncertainty in Table III includes: 1) uncertainty from the 
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TABLE III: Contributions to Experiment Uncertainty 



Contribution 


NS (mas/yr) 


WE (mas/yr) 


Statistical (modeled) 


16.8 


5.9 


Systematic (unmodeled) 






Parameter sensitivity 


7.1 


4.0 


Guide star motion 


0.1 


0.1 


Solar geodetic effect 


0.3 


0.6 


Telescope readout 


0.5 


0.5 


Other readout uncertainties 


< 1 


< 1 


Other classical torques 


< 0.3 


< 0.4 


Total Stat. + Sys. 


18.3 


7.2 



sensitivity analysis to the number of model parameters; 
and 2) small effects not incorporated in the model. 

We performed several important data analysis cross- 
checks. First, the gyro drift rate results of Table |LT| are 
confirmed by separate analyses of the segmented data. 
The 24 independent results from six data segments for 
each gyroscope are all consistent with the joint result 
within their confidence limits, demonstrating the inter- 
nal consistency of the model. Second, the misalignment 
torque coefficients k determined during the calibration 
phase proved to be in excellent agreement with the end- 
of-mission values estimated by both algebraic and ge- 
ometric analysis methods. No less impressive was the 
agreement between the time history k(t) throughout the 
mission obtained by the two methods. As for the roll- 
polhode resonance torques, the gyro dynamics model of 
Eq. ([IJ predicts that during a resonance the gyroscope 
orientation axis approximately follows a Cornu spiral. In- 
deed, that is typically observed in orbit-by-orbit gyro ori- 
entations determined by both data analysis methods. 

Further cross-checks from the geometric method in- 
cluded relativity estimates for 2 of the 4 gyroscopes, in 
both NS and WE directions, in statistical agreement with 
the algebraic results, and an estimate of the gravitational 
deflection of light by the Sun. IM Pegasi's closest ap- 
proach to the Sun occurs on March 11th, with ecliptic 
latitude 22.1° as compared with the grazing incidence 
value 0.265° in classical light deflection tests. The max- 
imum deflection predicted by GR is 21.7 mas; the ob- 
served 21 ±7 mas serves as a neat structural confirmation 
of Gravity Probe B. 

Lunar laser ranging has reported a measurement of the 
de Sitter solar geodetic effect to 0.7% [15] . Analyses of 
laser ranging to the LAGEOS and LAGEOS II spacecraft 
report a 10% - 30% measurement of the frame-dragging 
effect, assuming the GR value for the geodetic preces- 
sion [161 117) . GP-B provides independent measurements 
of the geodetic and frame-dragging effects at an accuracy 
of 0.28% and 19%, respectively. 
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